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Abstract 

Elastic production of vector mesons 7* N — > V N is the pomeron-exchange dominated 
diffractive reaction with much potential of probing the BFKL pomeron. The BFKL pomeron 
can conveniently be described in terms of the dipole cross section which is a solution of the 
generalized BFKL equation. In this paper we discuss, how the energy and Q 2 dependence 
of elastic production of vector mesons at HERA will allow scanning the dipole cross section 
as a function of dipole size r. We show that determinaton of the intercept of the BFKL 
pomeron requires measuring the p° and J/\l/ production at Q 2 ~ (100 — 200)GeV 2 and/or 
the quasireal photoproduction of the T. We present predictions for the effective intercept 
in the kinematic range of the forthcoming HERA experiments, which can shed much light 
on the nonperturbative component of the pomeron. 
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Elastic real and virtual photoproduction of vector mesons 

7 *p^Vp, V = p°, uj°, 0°, J/ty, T (1) 

is the typical diffractive reaction, dominated by pomeron exchange. Determination of pa- 
rameters of the BFKL pomeron [1], in particular of its intercept An>, is an outstanding 
problem for the HERA experiments on small x, and exploring the potential of reaction (1) 
is of great importance. In [2-4] we developed a novel approach to the BFKL pomeron in 
terms of the dipole cross section cr(£, r) which satisfies the generalized BFKL equation. Here 
r is the size of the color dipole, £ = \og[2.m p v / (Q 2 + my)] is the rapidity, my is the vector 
meson mass, Q 2 is the virtuality of the photon and v its energy in the proton rest frame. 
At very high v and/or very small Bjorken variable x one expects 

a{Z,r) = ajp{r)exp{AjpS) (2) 

and 

F 2 (x,Q 2 )ocexp(A IP 0~ (^) A " (3) 
with the intercept Ajp which does not depend on Q 2 , the real photoproduction Q 2 = not 
excepted. Understanding the onset of the BFKL regime i s the most pressing issue, 

because experimentally a steep rise of F2(x,Q 2 ) at HERA [5] is accompanied by a much 
slower rise of the real photoabsorption cross section [6]. 

What makes reactions (1) exceptionally important is the scanning phenomenon discov- 
ered in [7-11], by which elastic production amplitude probes the dipole cross section at the 
scanning radius 

r s « . . (4) 

y/m* v + Q 2 

The scale parameter C is rather large , and in [11] we gave an estimate C ~ 6. Changing 
Q 2 and the mass of the produced vector meson, one can probe cr(£, r) and measure the 
effective intercept A e jy(£,r) = d log cr(£, r)/d£ in a very broad range of dipole size r, from 
the nonperturbative regime of r > If down to the perturbative regime of r If. The 
HERA experiments are already amassing the data on production of vector mesons, and the 
purpose of this communication is to present the detailed predictions for reactions (0) from 
the approach to the BFKL pomeron developed in [2-4,12-15]. 



In order to set up the framework, we briefly review properties of the dipole cross section. 
In [3,4,12,13] we studied solutions of the generalized BFKL equation in the realistic model 
with the finite correlation radius R c for perturbative gluons and the running QCD coupling 
as(r), which freezes at large distances as{r > Rf) = aig = 0.8. With the R c = 0.3 / as 
suggested by the lattice QCD studies, Aip = 0.4 [4,12]. The major findings are: 

Firstly, the onset of the asymptotic behaviour (§,0) was found to be extremely slow. 
Namely, the effective intercept A e jf(£,r) varies quite strongly in the range of r and £ rele- 
vant to the HERA experiments on DIS. The usually discussed structure function F2(x,Q 2 ) 
receives contributions from a broad range of r and, for this reason, its x dependence does 
not allow a realiable experimental determination of the intercept Ajp. 

Secondly, by the diffusion property of the BFKL kernel [1,4,12], at asymptotic energies 
the behaviour (0) takes over at all values of r, including r > R c , thus unifying the asymptotic 
energy dependence of the bare pomeron-exchange contribution to hadronic scattering, the 
real photoproduction and deep inelastic scattering. Furthermore, by the same diffusion 
property, the intercept Ajp is almost entirely controlled by the region of r ~ R c [4,12]. 

Thirdly, there happens to exist a "magic" size r& ~ \Rc, at which the precocious 
asymptotic behaviour A e //(£, r) ~ An> will persist over the whole range of x at HERA. In 
[4,15] we showed how one can zoom at r ~ r&, measuring at HERA the charm structure 
function at Q 2 < 10 GeV 2 , the longitudinal structure function Fl(x,Q 2 ) at Q 2 ~ (10 — 
40)GeV 2 and the scaling violations dF T (x, Q 2 )/d\ogQ 2 at Q 2 ~ (2 - 10)GeV 2 . 

Fourthly, the BFKL pomeron describes the exchange by perturbative gluons with the 
typical interaction radius ~ R c . At low energy and large dipole size r 3> R c , the nonper- 
turbative scattering mechanism can take over. The rise of the BFKL perturbative cross 
section is driven by the rising multiplicity of perturbative gluons in the high-energy photon 
(hadron), for the nonperturbative mechanism one can expect the (approximately) energy- 
independent contribution a^ npt \r) to the dipole cross section. A large contribution from 
the <j( npt )(r) reconciles the steep rise of F 2 (x,Q 2 ) [5] with a much weaker energy depen- 
dence of the real photoabsorption cross section [6] . As we shall discuss below, the scanning 
phenomenon allows to study in much detail a transition between the nonperturbative and 
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perturbative regimes, which makes reactions (1) particularly important. 

The origin of the scanning phenomenon is as follows: In the dipole-cross section repre- 
sentation, the amplitude of the forward production 7* p — > V N of the (T) transverse and 
(L) longitudinally polarized vector mesons reads [7-11,16,17] 



2 _ N C C \ ■ v ' bo, 



J d 2 ra{t,r)j^ - {m\K Q (er)<j>(r, z) - [z 2 + (1 - zf^K^d^r)} 

dr 2 a(^r 



(m 2 v + Qiy 



/^^MQV) (5) 



, 2 . N c C vy /^E^2 y /Q s 



(2tt) 2 m v 



1 rdr 2 a(£,r 



/£^W,,*) (6) 



(m v + Q 2 ) 2 my 

Here N c = 3 is the number of colours, CV = |, | for the p°, u; , 0°, produc- 

tion, respectively, K 0jl (x) is the modified Bessel function, 

e 2 = m 2 q + z(l-z)Q 2 , (7) 

r) is the lightcone wave function of the qq Fock state of the vector meson in the mixed 
(z, ^-representation, where r is the transverse separation of the quark and antiquark and 
z is a fraction of the lightcone momentum carried by the quark. The normalization is such 
that 

^ / d 2 r [ -*?L- { m 2 2 ( r , z) + [z 2 + (1 - ,) 2 R0(r, z)] 2 } = 1 . (8) 

The dipole cross section can be related to the gluon structure function of the target proton 
G(x,q 2 ) [2,17,18] 



7T 2 A 

CT (£> r ) = —r 2 a s (r)G(x = exp(f ), q 2 = -j) oc r 2 
3 r z 



7-1 

exp(A IP , (9) 



_a s (r) 

where 7 = 4/3 Aip [14] and A 10 (the emergence of this large numerical factor is explained 
in [13]). The most important property of the dipole cross section is the color transparency 
driven dependence oc r 2 at small r. Because K 0jl (x) oc exp(— x) at large x, and the wave 



function of the vector meson is smooth, the amplitudes (|5||]) will be dominated by the 
contribution from r « 3/e, which in the nonrelativistic approximation of my ~ 2m g and 
2 ~ | leads to the scanning radius (f|) and the estimate C ~ 6. When the scanning radius is 
small, rs < Rv, where Rv is the vector meson radius, the amplitudes can be evaluated 

as 



1 

{m 2 v ~Q 2 ) 2 



M T oc rf <r(£, r s ) oc , 2 ; ^ , (10) 



m l oc r s a(£,r s oc oc M T , 11 

my my (m y + Q 2 ) 2 m v 

which must be contrasted to prediction of the vector dominance model oc l/(m v + Q 2 ). 
This departure from the vector dominance comes entirely from colour transparency property 
of the dipole cross section [11]. Production of the longitudinally polarized vector mesons 
dominates at Q 2 3> my. 

The scanning property is quantified by weight function Wt,l(Q 2 , r<2 ) m Eqs. (|5]||), which 
we introduced following [13]. The Wt,l(Q 2 \ r2 ) are sharply peaked (Fig. 1) and, because 
a(^,r)/r 2 is a slow function of r, see Eq. the position of the peak gives the scanning 
radius rs- At sufficiently large Q 2 such that rs < Rv, which typically means Q 2 ^> my, the 
position of the peak in the variable y = r 2 (Q 2 + m v ) very weakly depends on Q 2 . For pro- 
duction of the longitudinally polarized vector mesons, the peak corresponds to Cl(J/^) ~ 7 
and Ci{p°) ~ 8. For the transversely polarized p° mesons, with rising Q 2 the peak is slowly 
drifting towards Ct ~ 12 at Q 2 ~ lOOGeV 2 , which can be understood as follows [11]: At 
a fixed z, the scanning radius rs ~ 3/e. For the longitudinal photons, the wave function 
is peaked at z ~ ^, hence the estimate C ~ 6 is rather accurate (for the related discussion 
of the longitudinal structure function see [2,17,20]). In the opposite to that, the transverse 
structure function receives significant contribution from the very asymmetric pairs with 
z, (1 — z) ~ m 2 /Q 2 , which have a large, hadronic transverse size r ~ l/m q [2,17]. Because 
of the contribution of such asymmetric pairs, the effective scanning radius for the trans- 
versely polarized vector mesons decreases with Q 2 less rapidly than given by the simple 
estimate (^). The results shown in Fig. 1, were obtained with the relativized wave function 
of vector mesons, which has the correct short distance behaviour driven by hard gluon ex- 
change [16]. Notice, that because of very large numerical value of Ct,l, the scanning radius 



rs Eq. dJ) remains large, and the relativistic effects in wave functions are still marginal, in 
a very broad range of Q 2 . For instance, for the J/\& we find Ct ~ Cx and very close to the 
nonrelativistic estimate Ct,l ~ 6 even at Q 2 ^ rriy. 

At small Q 2 , the scanning radius rs is large, and one probes the dipole cross section in 
the nonperturbative domain of r. In Fig. 2 we present the decomposition of the dipole cross 
section into the nonperturbative component a^ npt \r) and the a^ pt \^^r) which is a solution 
of the perturbative BFKL equation^ : 

a{^r) = ^ npt \r) + ^ pt \^r) (12) 

This additivity of the bare pomeron cross sections and/or of the eikonal functions when the 
unitarization is considered, is the simplest assumption, and more refined treatment will be 
necessary in an analysis of the future high-precision data. The shape of the nonperturbative 
cross section is the largest unknown in the problem, and is mostly driven by our analysis of 
the Q 2 dependence of the photoproduction of vector mesons [8-11] and of structure functions 
at moderate Q 2 [18], which are well described by the dipole cross section calculated in [17]. 
Therefore, the tuned so as to reproduce the dipole cross section of ref. [17] at 

£ = £o — — log 0.03, which we take as the starting point for the BFKL evolution. In Fig. 2 
we show how the cr^(£, r) evolves with the rapidity £. At small r, the total dipole cross 
section is dominated by the BFKL cross section cr^(£, r). In this region, because of the 
so-called double logarithmic effects, the effective intercept A e jj(£, r) > Aip [3,4]. At large r, 
A e ff(l;,r) <C Aip because of the large contribution from the nonperturbative cross section 
0-( np *)(r). The analysis [4,13] has shown that there exists a magic radius r ~ r&, at which 
A e ff(£,r) ~ Aip starting already at moderate £. At this value of r, the nonperturbative 
cross section is small, and zooming at r ~ ta shall allow direct measurement of Ajp at 
HERA. 

The above discussion refers to the bare pomeron cross section. With the conventional 

Gaussian parametrization of the elastic scattering peak, da e i/dt oc exp(—B\t\), where \t\ is 
1 This Fig. 2 updates Fig. 1 of Ref. [14], in which the mistake was made in plotting cr( pt ) (£ , r) . This 
plotting mistake does not affect any of predictions made in [14], though. Our choice of cr^" p *- ) (r) in this 
paper is slightly different from that in [14]. 
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the momentum transfer squared, the profile function of elastic scattering takes the form 

™ = si (13) 

where b is the impact parameter. The rise of the bare pomeron cross section will eventu- 
ally conflict the s-channel unitarity bound T(b) < 1. Unique solution of the unitarization 
problem is lacking; for the crude estimate we apply the /C-matrix unitraization 

When the rising strength of the bare pomeron interaction T (b) ^> 1, the unitraized profile 
function tends to the black disc limit F(b) — > 1. The /C-matrix unitarization leads to a 
particularly simple form of the unitarized total cross section [19] 

^)(£,r) = 47r31og(l + ^) . (15) 

For a crude evaluation of the unitarization effects, we take the energy-independent diffraction 
slope B = 10GeV~ 2 . Evidently, the unitarization effects only become important when 
cr (£> r ) ^ 47r5 ~ 50mb, i.e., either at large r and/or at very high energy. As we shall see 
below, in the kinematical range of the HERA experiments, unitarization effects are marginal. 

Our predictions for the energy dependence of different photoproduction observables are 
shown in Fig. 3,4. In Fig. 3 we show the real photoabsorption cross section. It was calcu- 
lated using the wave function of the photon derived in [17], assuming the effective quark 
mass m Uy( i = 0.15GeV, which is the sole parameter in the quark wave function. Only the 
nonperturbative contribution to <Jtot{lP) is sensitive to this parameter. For the comparison 
purposes, we also show our predictions for <J to t{p°p), which is close to <r to i(7riV) and has a 
rise consistent with the observed trend of the hadronic total cross sections (for a recent 
review and high-energy extrapolations see [22]). We attribute the slow rise of (Jtot{lP) and 
<Jtot{p°p) to the large contribution from the nonperturbative large-size cross section (for the 
early discussion of such a scenario see [23]). The effect of the unitarization slowly rises with 
energy. It predominantly affects the nonperturbative cross section and, given the uncer- 
tanty in the absolute normalization of the cr( npt \r), at small energies it can be compensated 
for by upwards renormalization of the input a^ npt \r). The energy dependent part of the 
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unitarization correction is, however, the genuine effect and brings our results for (Jtoti.lv) to 
a better agreement with the experiment [6]. The simple /C-matrix unitarization, as well as 
the eikonal unitarization, lacks the so-called triple-pomeron contribution, which enhances 
the effect of unitarization at high energies ([2,19] and references therein) and can tame the 
somewhat too rapid a growth of cr tot (7p), a description of the low energy data [21] also can 
be augmented, by adding to the nonperturbative cross-section the Regge-behaving oc \ 
terms [22]. However, the purpose of the present analysis is understanding the gross features 
of the BFKL phenomenology rather than fitting the low-Q 2 experimental data. Recently, 
there was much discussion on the perturbative QCD contribution to the real photoabsorp- 
tion cross section (JtotilP), in which the eikonal is evaluated with the BFKL cross section 
substituted for the inclusive minijet cross section ([24] and references therein) which is 
incorrect (for the related criticism see [25]). 

In Fig. 4 we show how the unitarization affects the energy dependence of da/dt\ t =o for 
the real photoproduction of the p° and J/\&. In the former case the scanning radius is large, 
r s > If, in the latter case r s ~ 0.4f. We also show the unitarization effect for the virtual 
photoproduction of the longitudinally polarized p° meson at Q 2 = 120GeV 2 , appropriate for 
the scanning radius rs ~ ta ~ 0.1 5f. Notice a rapid decrease of the unitarization correction 
with the decrease of the scanning radius rs- 

For a crude estimate of elastic production cross section we can take B(p°) m lOGeV 2 
[6] and B(J/^) m 4GeV 2 [26]. For the real photoproduction of the p° at y/s = y/'2m p u = 
200GeV, we find a(ip — ► p°p) / &tot{lP) = 0.085 (for the bare, non-unitarized, cross sections 
this ratio equals 0.135 ), in good agreement with the ZEUS determination 0.1 ±0.04 [6]. For 
the real photoproduction of the J/^f at v — 150GeV we find a(jp — > J/^/p) ~ 16nb, which 
agrees with the E687 result of 17.9±4.0nb at v = 177GeV and 9.8±2.9nb at v = 121GeV [26]. 
The rise of the J/\I/ production cross section with energy observed by E687 is also consistent 
with our prediction in Fig. 4. Our estimate for the real photoproduction of the J/*f> at 
HERA at y/s = 200GeV, i.e., f = 8.3 is cj to t{lV -> J/^p) ~ 90nb, i.e., we predict the five- 
fold increase from FNAL to HERA. Theoretical calculations of the forward production cross 
section are more straightforward and are free of uncertainties with the energy dependence of 
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the diffraction slope, and it is very much desirable that the experimental data are presented 
in the form of da/dt\ t= o- More detailed comparison with the experiment will be presented 
elsewhere. 

In Fig. 5 we present our predictions for the effective intercept 

2 1 d\og(dcr/dt)\t=o) 



2 «9£ 

for the transverse and longitudinal polarizations of the J/\I/ and p°. The rise of A e ff(£, Q 2 ) 
with energy and Q 2 is an interplay of the BFKL perturbative and the nonperturbative cross 
section. On the other hand, for the J/\l/ and the longitudinally polarized p°, at Q 2 ~ 
C 2 jr\ ~ (100 - 200) GeV 2 we predict the precocious BFKL behaviour A e// (£,Q 2 ) « Ajp. 
Such a behaviour of A e ff(£, Q 2 ) is a very definitive prediction of our approach. Because of a 
large value of Ct{p°), in the transverse p° production the magic scanning radius r$ ~ ta and 
A e ff(l;,Q 2 ) « Aip are not attainable even at Q 2 ~ 200GeV 2 . This prediction is difficult to 
test, though, because o T <C ctl- It is interesting that one has r$ ~ r& also for the quasireal 
photoproduction of the Y(IS'), which reaction can even be favoured for the higher flux of 
quasireal photons. For the p° and J/^ production at large Q 2 , the kinematical range of 
HERA corresponds to exp(— £) « x > 10 -5 • (Q 2 /GeV 2 ), in which region the unitarization 
effects can still be neglected. We emphasize that precisely the same approach with the same 
dipole cross section, gives a very good description of the HERA data on ^(x, Q 2 ) [15]. 

Within our lightcone formalism, the dipole cross section is a universal quantity, and all 
the dependence on the process is contained in wave functions of the initial and secondary 
particle. The point we wish to make is that once the different production processes 7*p — > 
Vip are studied at values Q 2 so arranged as to have the same scanning radius rs Eq. (f|), then 
the production cross sections will exhibit identical energy dependence. One must compare 
the cross sections at energy z/, corresponding to the same rapidity £. For instance, we predict 
identical energy (£) dependence of the real photoproduction of the J/\I/ and of the virtual 
photoproduction of the longitudinal and transverse p° at Q 2 w 20GeV 2 and w 30GeV 2 , 
respectively. This rescaling from one vector meson to another constitutes an important 
cross-check of the whole formalism. 

More predictions can be made, which are specific of the small gluon correlation radius 



R c = 0.3f. Namely for this reason, the perturbative BFKL cross section cn p *)(£, r) is con- 
centrated at smaller r than the a^ npt \r), and the onset of the dominance of the BFKL cross 
section is followed by a change of the Q 2 dependence of the forward production cross section. 
In Fig. 6 we present our predictions for the energy dependence of 



D T = (m 2 v + Q 



2 i /-|2\4 d(T r p 



(It 



t=0 ' D i-( m v + Q) ■q^-^T 



t=0 



Q 2 d<7( 7 £-V T ) 



The gradual rise of Dl t with Q 2 comes predominantly from the rising density of gluons, 
see Eq. ([]). Notice, that the ratio do^jdaT rises much slower than Q 2 /rriy and exhibits 
a substantial energy dependence at fixed Q 2 , which can be tested at HERA. We strongly 
advocate studying the Q 2 dependence at fixed rapidity £ rather than the fixed energy v. 

Eqs. (P-|TT|) show that the vector-meson production amplitude measures the dipole cross 
section <j{£,,rs) and the gluon structure function of the proton. We wish to emphasize that 
this measurement is exceedingly sensitive to the r$ — Q 2 relationship. Typically, one will 
probe the gluon structure function G(^,Q 2 ) at the factorization scale q 2 ~ r(Q 2 + m v ), 
where 

t t ,l ~ (16) 

In the range of Q 2 ~ (10 — 100)GeV 2 of the practical interest, we find tt,l(J/^) ~ 0.2, 
t l{p°) ~ 0.15 and tt(p°) ~ (0.07 — 0.1). The emergence of such a dramatically small 
rescaling coefficients tt,l is a consequence of colour transparency. It was overlooked in [27]. 

A brief commment on the production of the 2S radial excitations p' and is in or- 
der. The wave function of the 2S state has a node, because of which the corresponding 
weight function WT^iQ 2 ,^ 2 ) changes the sign at r ~ Ry- The resulting cancellations of 
the r > Ry and r < Ry contributions to the production amplitude (the node effect) lead 
to a strong suppression of the p' / p° and ^/'/(J/^/) production ratio [7-11]. With increasing 
energy and with the onset of the dominance of a^ pt \^, r), which is concentrated at smaller r, 
the node effect will decrease. For instance, in the real photoproduction the \E f/ /(J/ v I / ) ratio 
is expected to increase by the factor ~ 3 from the CERN/FNAL energies to the highest 
energies available at HERA [28]. The interesting possibility [28] is that for the radially 
excited light mesons (0', u', p'), in the real photoproduction the V'/V ratio may exhibit the 
anomalous, nonmonotonic energy dependence, when this ratio first decreases with energy, 
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and then starts increasing. In the opposite to this, at large Q 2 such that rs <C Rv, the node 
effect disappears, and we expect p' / p° ~ ^f'/(J/^f) ~ 1 over the whole energy range. 

Conclusions: 

The purpose of this paper was to examine how the dipole cross section cr(£,r) is scanned 
in the virtual photoproduction of vector mesons. The BFKL dipole cross section, comple- 
mented by the nonperturbative dipole cross section at large r, gives the unified description 
of the photoproduction processes in the whole range of v and Q 2 , including the real pho- 
toproduction. The same dipole cross section gives a good description [14] of the proton 
structure function measured at HERA. Changing Q 2 , one can probe the energy dependence 
of cr(£,r) in a broad range of radii from the nonperturbative region of r ~ If down to the 
perturbative domain of r <C If. One can zoom at the magic radius and determine at 
HERA the pomeron intercept Aip, measuring the cross section of elastic p°, J/ty production 
at at Q 2 ~ (100 — 200)GeV 2 and the quasireal photoproduction of the T(IS'). We argued 
that this determination of A^. is not affected by the unitarity corrections. We make a strong 
point that it is the dipole size r rather than the photon's virtuality Q 2 , which controls the 
energy dependence of diflractive amplitudes: the value of Q 2 needed to scan the dipole cross 
section at the magic radius r& strongly varies from the process to process, c.f. the discussion 
in [15]. We predict very specific variation with energy of the Q 2 dependence of elastic pro- 
duction cross section and of the L/T ratio, which derive from the small correlation radius 
for perturbative gluons R c ~ 0.3f. Elastic production is found to probe the gluon structure 
function of the proton at an anomalously small factorization scale. 
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Figure captions: 

Fig. 1 - Weight functions Wt,l(Q 2 , r 2 ) as a function of a variable y = r 2 (Q 2 + m v ). 

Fig. 2 - Decomposition of the dipole cross section into the BFKL perturbative cr^(£,r) 
and the nonperturbative a^ npt " 1 (r) components. The growth of the BFKL cross section 
with energy (rapidity) is shown. 

Fig. 3 - The predicted energy dependence of the (B) bare and (U) unitarized total p°p and 
real photoabsorption cross section. The data shown are from the HERA [6] and FNAL 
[21] experiments. 

Fig. 4 - The predicted energy (rapidity) dependence of the (B) bare and (U) unitarized 
differential cross section of forward real and virtual photoproduction of the p° and 

Fig. 5 - The energy (rapidity) and Q 2 dependence of the effective intercept A e ff(^,Q 2 ) for 
the forward production of the transverse and longitudinal p° and J/ty. 

Fig. 6 - Predictions for the Q 2 dependence of the differential cross section of forward pro- 
duction and the ratio of the longitudinal and transverse cross sections. Shown are the 
quantities D TL and R defined in the text. 
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